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Abstract 

The differential-reduction algorithm, which allows one to express generalized 
hypergeometric functions with parameters of arbitrary values in terms of such func- 
tions with parameters whose values differ from the original ones by integers, is 
discussed in the context of evaluating Feynman diagrams. Where this is possible, 
we compare our results with those obtained using standard techniques. It is shown 
that the criterion of reducibility of multiloop Feynman integrals can be reformu- 
lated in terms of the criterion of reducibility of hypergeometric functions. The 
relation between the numbers of master integrals obtained by differential reduction 
and integration by parts is discussed. 
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1 Introduction 



It is commonly accepted that any multiloop and/or multileg Feynman diagram in 
covariant gauge within dimensional regularization [1] may be treated as a generalized 
hypergeometric function-*^ [3]. Starting from its a representation, any Feynman 
diagram may be written in the form of a Mellin-Barnes integral [4,5], 



where Ya are algebraic functions of external kinematic invariants and A,B,C,D are 
some matrices depending in a linear way on the dimension of space-time n, which 

is an arbitrary complex number, and the powers of the propagators. 

By the application of Cauchy's theorem,^ this integral can be rewritten as a 
linear combination of multiple series,^ 

For real diagrams, some of the variables may be complex number. We call this 
type of variable a "hidden variable" and the corresponding index of summation a 
"hidden index of summation." In all existing examples, the representation of Eq. (2) 
belongs to a Horn-type hypergeometric series [8] if the hidden index of summation 
is considered as an independent variable. 

For the reader's convenience, we recall that the multiple series Y1a=o C{m)x^ is 
called Horn-type hypergeometric if, for each ? = 1, . . . , r, the ratio C(m + ei)/C{fn) 
is a rational function in the index of summation (mi, • • • , rur) [8,9]. The coefficients 
of such a series have the general form 

c{rn) = n Ar^(^) iii7V, !^ ' 

where N,M > 0, Xj,Sj,jj are arbitrary complex numbers, : ^ Z are 

arbitrary integer- valued linear maps, and R is an arbitrary rational function [8,10]. 

However, to our knowledge, the proof that any Feynman diagram can be de- 
scribed by a Horn- type series does not exist. There is another way to proof this 
statement. The classical a representation of a Feynman diagram is a particular case 
of the generalized Euler integral representation of the Gel'fand-Kapranov-Zelevinski 
system [8] that is related to the Horn-type series representation. We call the Feyn- 
man diagram representation of Eq. (1) or the equivalent representation of Eq. (2) a 
hypergeometric representation [8]. 

For Horn-type hypergeometric functions, there are so-called step-up (step-down) 
operators i^x) [9)11]- These are differential operators which, upon application 



^This statement is also valid for phase-space integrals (see, e.g., Rcf. [2]). 

^This is true if all arguments of the F functions are different. Otherwise, an additional regularization 
for each propagator (or the introduction of extra masses) is necessary. See, e.g., Ref. [6]. 

^One of the first examples of this type of representation for a Feynman diagram was given in Ref. [7]. 
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to a hypergeometric function Sa, shift the value of one its upper (lower) param- 
eters by unity, as HI^Sx = Sx+i {H^Sx = Sx-i)- Takayama [12] proposed an 
algorithm that allows one to construct inverse differential operators, the step-down 
(step-up) operators (-^a )' starting from direct operators and systems of dif- 
ferential equations for hypergeometric functions. These operators satisfy the re- 
lations B^Sx+i —5- Sx {B^Sx — 5" 5'a+i)- Takaj^ama pointed out [12] that inverse 
operators are uniquely defined for any hypergeometric function with an irreducible 
monodromy group, which implies that the parameters and the differences between 
upper and lower parameters are not integer. By the action of such differential op- 
erators on a hypergeometric function, the value of any parameter can be shifted 
by an arbitrary integer. We call this procedure of applying differential operators 
to shift the parameters by integers differential reduction.'^ An important step of 
Takayama's algorithm is the construction of a differential Grobner basis for the 
system of differential equations for hypergeometric functions.^ 

It is quite surprising that the technique for the reduction of Feynman diagrams 
advocated here, namely to split a given Feynman diagram into a linear combination 
of Horn- type hypergeometric functions with rational coefficients and to subsequently 
apply differential reduction, has never been elaborated and that its interrelation with 
the well-known intcgration-by-parts (IBP) technique [15] has never been discussed. 

The aim of this paper is to demonstrate how the differential-reduction algorithm 
may be successfully applied to evaluate Feynman diagrams. For simplicity, we con- 
sider here only the particular case of Horn-type multiple hypergeometric functions, 
i.e. the functions p^iFp, and some Feynman diagrams with arbitrary powers of 
propagators, which are expressible in terms of these functions. 

The structure of this paper is as follows. In Section 2, we discuss the differential- 
reduction algorithm for generalized hypergeometric functions pj^iFp. This is a cru- 
cial step towards proving theorems on the construction of all-order e expansions 
presented in Refs. [16-19]. Section 3 illustrates the application of differential re- 
duction to several Feynman diagrams of phenomenological interest. In Section 4, 
we discuss how the counting of master integrals in the differential-reduction ap- 
proach is related to that in the IBP technique. The results of our analysis are 
briefly summarized in Section 5. In Appendix A, we describe interrelations between 
a set of basis functions generated by the differential-reduction algorithm and a set 
of hypergeometric functions whose higher-order e expansions were constructed in 
Ref. [20]. 

'*Thcre are various publications on contiguous relations for hypergeometric functions, i.e. algebraic 
relations between hypergeometric functions of several variables with shifted values of parameters, starting 
from the classical paper by Gauss [13]. To our knowledge, a closed algorithm for the algebraic reduction 
of Horn-type hypergeometric functions of several variables, i.e. an algorithm for solving these algebraic 
relations, does not exist. 

^The idea of using the differential-Grobner-basis technique directly for the reduction of off-shell Feyn- 
man diagrams, without spitting them into linear combinations of hypergeometric functions, was proposed 
by Tarasov [14]. 
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2 Differential-reduction algorithm for the gen- 
eralized hypergeometric function p+ii^p 

2.1 Notation 

Let us consider the generalized hypergeometric function pFq{a;b; z), defined by a 
series about 2; = as 



pFq{d;b;z) = pFq 



nLi(«i)fe 



(4) 



where {a)k = r(a + k)/T{a) is the Pochhammer symbol. The sets a = (ai, • • • , Op) 
and h= (61, • • • ,hq) are called the upper and lower parameters of the hypergeometric 
function, respectively. In terms of the differential operator 0, 



dz ' 



(5) 



the differential equation for the hypergeometric function pFq can be written as 



z\{{e + ai)-e\{{e + hi-i) 



1=1 



pFq{d;b;z) = 0. 



(6) 



Hypergeometric functions which differ by ±1 in the value of one of their parameters 
are called contiguous, and linear relations between contiguous functions and their 
derivatives are called contiguous relations.^ In particular, the following differential 
identities between contiguous functions are universal [11]: 

pFq{ai + 1, a; b; z) = B^pFq{ai,a; b;z) = — {6 + ai) pFq{ai,d; b; z) , 

pFq{a;bi - l,b;z) = H^^pFq{a;bi,b; z) = _ {6 + bi- 1) pFq{d;bi,b; z) , (7) 

which directly follow from the series representation of Eq. (4). The operators 5+ 
(i^^ ) are the step-up (step-down) operators for upper (lower) parameters of hyper- 
geometric functions. 



2.2 Non-exceptional values of parameters 

In Ref. [12], it was shown that, for given step- up (step-down) operators, inverse 
step-down (step-up) operators, which are uniquely defined modulo Eq. (6), can be 
constructed. This type of operators were explicitly constructed for the hypergeo- 
metric function p+iFp by Takayama in Ref. [22]. For completeness, we reproduce 

full set of contiguous relations for an arbitrary hypergeometric function pFg was considered by 
Rainville [21]. For pFg, there are 2{p + q) contiguous functions. 
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his result here: 

p+iFp{ai - 1, a; b; z) = B^^p+iFp{ai,a; b; z) , 



b: 



04 
Ci 



ti(x) 



-aiY[{bj - l- tti) , 

xlfj=i{x + bj -1)-Ci 

X + ai 



(8) 



p+iFp{d- bi + 1, b; z) = H+p+iFp{d; bi,b; z) , 

h-l 



Y[{9 + bj-l)-sm 



p+1 

di = + aj - bi) , 

miix + aj)-di 



x + bi- 1 



(9) 



where \a^a+i means substitution of a by a + 1. The functions ti(x) and Si(x) are 
polynomials in x, so that Eqs. (8) and (9) are polynomials in the derivative 6. Let 
us introduce the symmetric polynomial Pl^\{rk}) as follows: 



p p 



j=0 j=0 



(10) 



so that 



k=l 



p^\{rk})= t: n ^ 

ii,--- ,v=l ii<---<ij 



(11) 



For example, we have Pi'\{rk}) = Z]^=i ''j £^nd Pp'\{rk}) = 11^=1 r^- Then, we 
may write 

P o+l f P j 



tiix) = Y.P^%{{br-l}) 



j=0 



X + Qi 



Y.Pp-j(ibr - l})Y.x^-H-ai)' . (12) 

j=0 k=0 



A similar consideration is valid also for the last relation in Eq. (9): 



j=0 



fc=0 
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The differential reduction has the form of a product of several differential step-up 
and step-down operators, and , respectively:^ 



F{a + m;h + n; z) 



a}) 



±J F{a-Mz), 



(14) 



so that the maximal power of 9 in this expression is equal to r = rrii + rij . 
In a symbolic form, this may be written as 

i?(aj, bj,z)F{a + m;b + n;z) = [51(02, bj,z)9^ + •••-!- Sr+i{ai, bj,z)] F{a; b; z) ,(15) 

where R and {S arc some polynomials. Since the hypergeometric function p+ii^p (a; 6; z) 
satisfies the following differential equation of order p+1 [see Eq. (6)] : 



{l-z)9P+\+^Fp{a-b-z) 



.r-=l 



p+1 



= E - Plti-m - 1})] 0^ + zY{aA p+,Fp{a; b; z) , (16) 



fe=i 



it is possible to express all terms containing higher powers of the operator 6, 6^ with 

k > p+1, as a linear combination of rational functions of z depending parametrically 
on a and h multiplied by lower powers of 9, 9^ with j < p. In this way, any function 
p+i-Fp(a + m]b + k] z), where fh and k are sets of integers, is expressible in terms of 
the basic function and its first p derivatives as 

S{ai,bj,z)p+iFp{a + m-,b + k-,z) = {Ri{ai,bj, z)9^ + R2{ai,bj, z)9^''^ 
+ \- Rp{ai,bj, z)9 + Rp+i{ai,bj, z)} p+iFp{a;b; z) , (17) 

where 5" and Ri are polynomials in the parameters {ai} and {bj} and the argument 
z. 

From Eq. (8) it follows that, if one of the upper parameters aj is equal to unity, 
then the application of the step-down operator to the hypergeometric function 

p+iFp produces unity, p^iFp{l, a; b; z) = 1. Taking into account the explicit form 
of the step-down operator B^ ^ 



B 



1 



WUih - 1) 



JJ(5^._l) + ^pW^.({5,-l})0^- 



, (18) 



we obtain the differential identity 
( 



> p+iFp{l,d; 6; z) 



+ < 



i:[pi%i{h-l})-zPi%{{a,}) 



}p+iFp{l,a;b;z) = l[{bk - 1) .(19) 
k=l 



'^Due to the relation pFq{a; b; z) — 2: 11^=1 '^j/ (111=1 ^j)p-^«('^ + 1;^+ ^t^), not all step-up oper- 
ators are independent. In fact, {UUK^(Ult\B+JF{a;b;z) = iUltlB+JiUUK^^^^'^^'^) = 
nU^i/{mt\ak){d/dz)F{d-,b-,z). 
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The case where two or more upper parameters are equal to unity, e.g. ai = 02 = 1, 

does not generate any new identities. As a consequence, if there is a subset I of 
positive integers in the set of upper parameters, the reduction procedure has the 
modified form: 

P{ai,bj,z)p+iFp{l,a + rh-,b + k;z) = Ri{ai,bj,z) + ^R2{ai,bj, z)0^~^ 
+ \- Rp{ai,bj,z)9 + Rp+i{ai,bj,z)^ p+iFp{l,d-b;z) . (20) 

Let us write exphcit expressions for the inverse operators for several hypergeo- 
metric functions. For the Gauss hypergeometric function 2F1, we have: 



2F1 



ai — 1, a2 
bi 



1 



bi - ai 



[(1 - z)e + 61 - ai - 022;] 2F1 



ai,a2 
bi 



(21) 



2F1 



ai,a2 
bi + 1 



bi \ , d , 

[bi — ai){bi — 02) [ dz 



2F1 



ai,a2 
bi 



For the hypergeometric function 3F2, the inverse differential operators read: 

= {(1 - z)e^ + [(61 + 62 - 1 - ai) - z{a2 + as)] 9 + (bi - ai)(62 - oi) - 2:0203} 
X 3F2 



01,02,03 
bi,b2 



3-^2 



bi 



01,02,03 
bi + 1, 62 

I — Z „r, 



z] {ai- 6i)(o2 - 6i)(o3 - 61) 
62-1 

(oi + 02 + 03 - 61) 



— [010203 - (01 - 6i)(o2 - 6i)(o3 - 61)] } 3F2 



01,02,03 
bi,b2 



(22) 
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For the hypergeometric function 4F3, the differential operators read: 
ai — 1,02,03,04 



z\{bi- oi)(62 - oi)(63 - 01) 



= {(1 - z)e^ + [(61 + 62 + 63 - 2 - Ol) - Z(02 + 03 + 04)] i 



+ 



+ &1&3 + &2&3 + of + (Oi + 1) 1 - ^ 6j - z(0203 + O2O4 + O3O4) 



4i^3 



+(61 - oi)(62 - ai){h2. - oi) - 2:02030414^3 

«!, «2, «3) «4 



«i,a2,a3,a4 



61 



+ 



hi + 1,62,^3 
' 1 - z 



+ 



zj (01 - 6i)(o2 - 6i)(o3 - 5i)(a4 - 61) 
"62 + &3 -2 



(oi + 02 + 03 + 04 — 61) 



(b2-l)ib3-l) , ^ (yr- 
(01O2 + 01O3 + 0104 + 02O3 + 0204 + 0304) + bi Qj — 

^ \j=i 
^ [01020304 - (oi - 6i)(o2 - 6i)(o3 - 6i)(o4 - 61)] I 4F3 ( "'^^'^"^ 



bi 



(23) 



2.3 Differential-reduction algorithms for special values 
of parameters 

A special consideration is necessary when 6, = Oj + 1. In this case, the expressions 
in Eqs. (8) and (9) are equal to zero. Let us start with the case when all Oj are 
different. In this case, the following relation should be applied (see Eq. (15) in 
Chapter 5 of Ref. [21], Chapter 7.2 in Ref. [23], and Ref. [24]): 

o,fe, ••• 



(o - b) pFq 



+ 1,6 + 1, 



bp-lFq^l 



o, ••• 

+ 1,' 



(24) 



Repeating this procedure several times, we are able to split any original hypergeo- 
metric function with several parameters having unit difference into a set of hyper- 
geometric functions with only one kind of parameters having unit difference (for the 
particular cases, see Eqs. (7.2.3.21)-(7.2.3.23) in Ref. [23]): 

^ f {ai}ri,{a2}r2r-- ){am}rr„,Ci,--- ,Cfe 



where {0^}^^ denotes Vj repetitions of Oj in the argument list, 

m 

R = '^rj , Cij^ai, 6j / 1 + Oj , Cj j^l + bj 



(25) 



(26) 
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For a special set of parameters, Eqs. (7.2.3.21) and (7.2.3.23) in Ref. [23] are useful: 



TP 
pJ^q 



{a}p-2,p,cr 
{b}q-2,P + n,a + 1 



z = 



(p-(T- 1)^ 



p-lFq-1 



{P)n p 
[P - 0-)n 

{a}p-2,P 
{h}q-2,P + k 



{a}p_2,CT 

{6}<,_2,(T+1 



(27) 



(m,-l)!^=i ^1=0 ■■■^"=0 

7 = 1 J 



n 

1=1 



n 



o-j + ii - o"ifc - iik 



p—n+lPq—n+1 



{a}p_„,(Tjk + jfc 

{6}5_„,(7fe + jfe + 1 



• (28) 



In Eq. (28), m„ are integers, all Ui are different, and, ii Ui — = N with A?^ 
1, 2, • • • , then < AT. 

Let us return to the last expression in Eq. (25) and rewrite it as 



^ {a + m}r, a-i + ki, - ■ ■ , ap+i_j. + /cp+i- 
P+i^P \^ {1 + a + m}q, 61 + /i, • • • , + /p-, 



(29) 



where m, {A;^}, and {Zj} are integers and a, {a^}, and are parameters of the 
basis function. Using the reduction procedure described in Sec. 2.2, we may convert 
this function as 



P+^^P\{l + a + m}q,{hu+'lk]p-q 



{a + m}r, {aj + kj} 



p+l— r 



^ {a + m}r, {aj + m}p+i_r 
P+i^P \^ {1 + a + m}q, {bk + m}p_g 



z , (30) 



and then apply the differential relation (see Eq. (7.2.3.47) in Ref. [23])^ 



F 

pJ^q 



{ai + m]p 
{hk + m}g 



ULdihU} ( d 



— r.F, 



p-L q 



{ai}p 

{h}q 



(31) 



where the derivative d/ (dz) could be rewritten in terms of 6, with the help of 

1' 



d 

dz 



e{e -!)■■■ {e -m + 1) . 



(32) 



If m > p, Eq. (31) can be converted to a differential identity of order p—1. Putting 
everything together, we obtain 



P(n, h- y\ .F I {a + ^}r,{aj + kj}p+l- 

= mak,hj,z)e^ + ■■■ + Qp+,{au,h^,z)U,Fp (^^^ + ^aCK ■ ■ ■ X-.-i ') ' ^^^^ 



^Another useful relation is Eq. (7.2.3.50) in Ref. [23]: 

[O'ilp 



.(7+n-l p ( {ai}p-i,o- 

P Q 



{ai}p-i,a + n 
{h}q 
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We remark that Eq. (33) can be further simphfied by using Eq. (7), which can be 
written as 



A, a 
l + A,b 



fc-1 



j=0 

Recursive appHcation of this expression allows us to reduce higher powers of deriva- 
tives in Eq. (33) as 



z] . (34) 



A 



{A}r,a 

{1 + A}r,b 



j=0 



p+i-jFp-j 



{l + A}r-j,h 



z] , (35) 



where q < r. 

For a particular set of parameters, a further simplification of Eq. (33) can be 
achieved. For example, for a = 1, we have 



p-Tq 



1) {0'i}p-l 

2,{6fc}g-l 



p-lFq-l 



{tti - l}p-l 



z -1 



, (36) 



where aj,bk ^ 1. 

Further useful relations for particular values of parameters of the basis function 
were derived in Ref. [24], namely 

{a}p,b 



p+i-Fp 



{1 + a}p 



V(jp) Jo {1-zty 



and, for 6 = 1, 



p+i-t^P 



1, a, • • • , a 
a + 1, • • • , a + 1 



z] =aP^{z,p,a) , 



oo u 



where ^(z,^, a) is the Lerch function defined as [9] 

oo 

$(z,p,a) = 

so that 



1 /"°° e~(""^)* 



k=0 



(a + k)P T{p) Jq e* - z 



Li„ (z) = z^{z,n, 1) . 



(37) 
(38) 

(39) 
(40) 



2.4 Criteria of reducibility of hypergeometric functions 



In this section, we formulate the criteria of reducibility of the hypergeometric func- 
tion pFq{a;b;z), i.e. we state under which conditions the hypergeometric function 
pFq{a; b; z) and its derivatives are expressible in terms of hypergeometric functions 
of lower order and/or with lower derivatives. 

We call the result derived by Karlsson [25] the first criterion of reducibility of 
the hypergeometric function pFq{a; b; z) to hypergeometric functions of lower order. 





+ mi, ■■■ ,bn + m„,a„+i, • 


\ 


) bn, bn+1, ■ ■ ■ ,bq 


mi 


m„ 




• ^ A{n,---jn)z-^"p-nF^ 


ji=0 


jn=0 



q—n 



bn+1 + Jnr 



dp ~\~ Jn 
, bq + Jji 



(41) 
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where rrij are positive integers, Jn = ji + ■ ■ ■ + jm and 



mi 
ji 



rUn 
3n 



(62 + ^2) Ji (63 + "is) J2 ■■■{K + rnn)j„_i (On+l) j„ • • • (flp) J„ 



(bi)ji (62)72 ■ ■ ■ {K)jAK+i)jn ■ ■ ■ {bq)j„ 



(42) 



In particular, we have 
/61 +mi,a2,-- - ,aj 

61,62, •••6, 



mi 

Y. 

j=0 



, /mi \ (02) 



\ j J (6i)i • • • {bq)j 



02 + i, • • • Op + j 

62 + j, • • • 6g + J 



In words, the first criterion of reducibihty of the hypergeometric function reads: 

Criterion I 

The hypergeometric function pFq(a;b;z) which has pairs of parameters 
satisfying ai = bi + rrii with rrii being positive integers is expressible in 
terms of functions of lower order according to Eq. (41). 

Equations (24) and (25) yield the second criterion of rcducibility of the hyperge- 
ometric function pFq{a; b; z) to functions of lower order. In its explicit form, it was 
derived in Ref. [24] (see Eqs. (18)-(20) in Ref. [24]). Assuming that (i) {ai, • • • a„} 
arc different and (ii) if a,. — = N with N = 1, 2, . . ., then rrii < we have (see 
also Eqs. (7.2.3.21) and (7.2.3.23) in Ref. [23]) 



^ ^ Vai + l + mi,--- ,a„ + 1 + m„,6„+i, • • • ,6^ 



n mi 

\ ^ \ ^ {-mi)j 



n 



1 




3) 



mr + 1 



\ai + \+ j,hn+i,--- ,6, 



(44) 



so that we may formulate it as: 
Criterion II 

The hypergeometric function pFq{a; b; z) which has two or more pairs of 
parameters satisfying &j = + mj -|- 1 with rrii being positive integers is 
expressible in terms of functions of lower order if additional conditions 
on the parameters ai (see Eq. (44)) are satisfied. 

Equations (34) and (35) are considered as the third criterion of reducibihty: 

Criterion III 

The result of the differential reduction of a hypergeometric function of 
the type pj^iFp{A + rh, a + k; 1 + A + m, b+l; z), where m, k, and I are sets 
of integers, is expressible in terms of the function p^iFp{A, a; 1 + A, 6; z) 
and hypergeometric functions of lower order and their derivatives. 

Equation (19) constitutes the fourth criterion of reducibihty: 
Criterion IV 

If one of the upper parameters of a hypergeometric function is an integer, 



(43) 
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the result of the differential reduction of this hypergeometric function 
has one less derivative and is described by Eq. (20). 

Criteria I-IV of reducibility of hypergeometric functions are mucii simpler ttian 
the ordinary criterion of reducibility of Feynman diagrams [26] . 



2.5 All-order e expansions of hypergeometric functions 

Recently, several theorems on the structure of the coefficients of all-order e ex- 
pansions of hypergeometric functions about integer and/or rational values of their 
parameters have been proven [16-19,27-30]. For a recent review, see Rcf. [30]. 

In this paper, we mainly consider hypergeometric functions of the type of Eq. (109). 
There still does not exist a rigorous mathematical proof regarding the structure of 
the coefficients of the all-order e expansions of hypergeometric functions of this 
type beyond the Gauss hypergeometric functions 2F1 [16,19]. Through functions of 
weight 4, the coefficients of the e expansions were constructed in Ref. [20,31], as 
discussed in details in Appendix A), which is sufficient for next-to-next-to-leading- 
order (two-loop) calculations. Recently, more coefficients of the £ expansions of the 
Clausen functions 3F2 have been evaluated in Ref. [32]. 



2.6 Reduction at z = 1 and construction of £ expansion 

The value 2: = 1 is a particular case of a "hidden" variable. It is evident that 
the application of Eq. (6) to the r.h.s. of Eq. (17) gives rise to the generation of 
factors 1/(1 — z)'^, so that the direct limit z — > 1 cannot be taken. Let us recall 
that the hypergeometric series defined by Eq. (4) converges [9] for \z\ = 1, when 
Rc^^^]^ bj — RcJ2^^l aj > 0, so that, if the l.h.s. of Eq. (17) converges, the r.h.s. 
of the equation does also exist. 

The main idea is to convert the original hypergeometric function p_|_iFp(a; b; z) to 
a function of argument \ — z. However, beyond type 2F1, the analytical continuation 
of the hypergeometric function pj^iFp{A\B]z) is not expressible in terms of 

functions of the same type, but has a more complicated structure [33]. Nevertheless, 
we can perform the analytical continuation z — ^ 1 — z of the coefficients of the e 
expansions of hypergeometric functions entering the r.h.s. of Eq. (17). It was shown 
in Ref. [34] that, under the transformation z 1— z, hyperlogarithms are expressible 
again in terms of hyperlogarithms. In this way, if the coefficients of the e expansion 
are expressible in terms of hyperlogarithms, there is an opportunity to find the limit 
z ^ 1 — z oi the differential reduction, but only in fixed orders of the e expansions 
and only for such values of parameters of the hypergeometric functions for which 
the analytical structures of the coefficients are known. 

An alternative approach to evaluate hypergeometric functions at z = 1 was 
discussed in Refs. [22,27]. 
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3 Application to Feynman diagrams 



As an illustration of the differential-reduction algorithm, let us consider the dia- 
grams shown in Figs. 1 and 2.^ In general, Feynman diagrams suffer from irreducible 
numerators. Using the Davydychev-Tarasov algorithm [37], any tensor integral may 
be represented in terms of scalar integrals with shifted space-time dimensions and 
arbitrary (positive) powers of propagators. In our analysis of the structures of the 
coefficients of the e expansions, we distinguish between two cases corresponding to 
even value n = 2m — 2£ and odd value n = 2m — 1 — 2£ of space-time dimension, 
where m is an integer. 











Figure 1: One- loop vertex diagrams expressible in terms of generalized hypergeometric 
functions. Bold and thin lines correspond to massive and massless propagators, respec- 
tively. 



Since all diagrams shown in Fig. 2 contain massless subloops, we present for 
completeness the result for the g-loop massless sunset-type propagator with q + 1 
massless lines. It is given by 



-^0(^^,0-1, 0-2, • • • ,(^q+l] 



[A;f [A;|]<^2 • • • [k^]"^ [{ki + k2 + --- + kg+ p)2]<^.+i 



'q+l 



r(§-^ik)l [zi-v/2]^r(a-fg) 



.ifc=i 



r(§(g + l)-a) 



(p2)f'/- ,(45) 



few examples of the differential-reduction algorithm applied to the reduction of Feynman diagrams 
were presented in Refs. [35,36]. 
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Figure 2: Diagrams considered in the paper. Bold and thin lines correspond to massive 
and massless propagators, respectively. 



where 

s 

(T = J2c7k, (46) 

fc=l 

with s = q + 1. Some of the massless lines can be dressed by insertions of massless 
chains. In this case, the respective ak can be represented as ak = Vk — Rk^, where 
rk and Rk and integers. 

3.1 One- loop vertex: particular cases 

We wish to remind the reader that any one-loop vertex diagram with arbitrary 

masses, external momenta, and powers of propagators can be reduced by recurrence 
relations, derived with the help of the integration-by-parts technique [15], to a ver- 
tex master integral plus propagator master integrals and bubble integrals (with all 
powers of propagators being equal to unity), or, in the case of zero Gram and/or 



14 



Cayley determinants, to a linear combination of propagator and bubble integrals. 
In terms of hypergeometric functions, vanishing Gram and/or Cayley determinants 
correspond to the situation where the hypergeometric function describing the orig- 
inal vertex diagram can be reduced to a Gauss hypergeometric function with the 
following sets of parameters [29,38,39]: 



2F1 



/2 + f 



2-^1 



/2 + I 



2 ' 



(47) 



where {la} are arbitrary integers. In the case of non-zero Gram and/or Cayley de- 
terminants, the one-loop vertex master integral is expressible as a linear combination 
of Gauss hypergeometric functions and Appell functions Fi [40]. 

For the diagrams shown in Fig. 1, both the Gram and Cayley determinants are 
non-zero. The hypergeometric representation of the corresponding master integral 
was published in Ref. [30]. 

• Diagram Ci with arbitrary powers of propagators is expressible in terms of 
two hypergeometric functions 3F2 as (see also Eq. (3.44) in Ref. [41]) 



Ci 



(-m2)§-^i23 



7r2 



r(ji23-f)r(§-ji2) 



r (!) m) 



-3^2 



J123 - 2) Jl! J2 
2>1 + .712 - 2 



X3F2 



n 

2-J12 



r(§-Ji)r(f-j2)r(ji2-t) 
r(n-ii2)r(ii)r(i2) 

J3,|-ji,f-i2 Q^' 

n-il2,f -jl2 + 1 




(48) 



Here and in the following, we use the short-hand notations jab = ja + jb and 
jabc = ja+jb+jc for compactness. In accordance with the differential-reduction 
algorithm, each of the two 3F2 functions in Eq. (48) is expressible in terms of 
a 2F1 function with one unit upper parameter, namely 



2F1 



1,1 



2F1 



l,Ii + f 
h + n 



(49) 



respectively, plus rational functions of z. Standard approaches yield one vertex 

master integral plus massive bubble and massless propagator integrals. The 
latter two types of integrals are expressible in terms of products of Gamma 
functions and correspond to the rational functions in our approach. 
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For diagram C2 with arbitrary powers of propagators, the result is 
C2 



(— m 



2w-jl23 



r (ii23-i) r (f -ji2) r {n-h^-2j2) 



+ - 



r(n-ii23)r(f-i2)r(i3) 

n 

Q2^2-^-r(§-ii)r(ji2-f)r(§-j23) 



J123— 2' Jl' J2 

l+ii2-f,l+ii3+2j2-n 



X3-^2 



r(n-ii23)r(ii)r02) 



J3, 2 -J2, 2 -Jl 

f -ii2 + i,i+i23-i 




(50) 



Similarly to the previous case, each of the two 3F2 functions in Eq. (50) is 
expressible in terms of a 2^1 function with one unit upper parameter, namely 



1,1 
I\ — n 



1,/i + i, 

/2- 2 



(51) 



respectively, plus rational functions. Standard approaches yield one vertex 
master integral plus massive bubble and massless propagator integrals. 

• Diagram C3 with arbitrary powers of propagators is expressible in terms of 
one 4F3 function as 



C3 



2^^-..3^(Jl23-i)^(f-J3) 



ii-"7r2 



= (-m^)2-^i23 



r(ji2)r(f) 



1F3 



/ J123 - f , jl, j2, f - js 
I n 212 ji2+l 

\ 2' 2 ' 2 



Am?' 



(52) 



{1,^1X3^2 



(53) 



In accordance with the differential-reduction algorithm, this function may be 
written in terms of a 3F2 function with one unit upper parameter and its first 
derivative, 

^1 - §,f - 1 + ^2 

§ +I2, 5 +/3 

and a rational function. Standard approaches yield one vertex and one prop- 
agator master integral plus massive bubble integrals. 

Diagram C4 with arbitrary powers of propagators is expressible in terms of 
one 3F2 function as 



C4 



10 

2\-9-J123 



r(ji23-i)r(f-ji) /^■,23-f,ii,j2 



r(i23)r(§) 



3-£^2 



2>J23 



^ ) • (54) 



zi-"7r2 

In accordance with the differential-reduction algorithm, this function may be 
written in terms of a 3F2 function with one unit upper parameter and its first 
derivative, 

7l + l,§+/2 

and a rational function. Standard approaches yield one vertex and one prop- 
agator master integral plus massive bubble integrals. 



{l,^}X3i^2 



(55) 
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3.2 Ef2o 

Let us consider the g-loop bubble diagram £^^20 depicted in Fig. 2 with two massive 
hnes of different masses and q — r plus r massless propagators, 

Ei2Q{rir?,M'^,a,P,ai, - ■ ■ , CTg-^, Pi, • • • ,Pr) = 



I 



1 



(56) 



The special case q = r was analyzed in Ref. [35]. The first non-trivial diagram of 
this type corresponds to a three-loop bubble {q = 3) with r = 1. The Mellin-Barnes 
representation of Eq. (56) reads: 



Ei2o{m^,M'^,a,P,ai,--- , C7q_r, pi, • • • ,Pr) 



• "I — 



r(a)r(/3)r (§) 



^ r(a,)' r(pfc) j r(f(g-r)-a) 



r(§-P.) l r(a-|(g-r-l)) /-^^^rn^ 



2 \ t 



r(-t)r (p-§(r-l)+t) r {a+p-^r+t) r (f (g-r)-(7+i) T (g+^g- §(g-r)-t) 
where cr is defined by Eq. (46) with s = q — r and 

r 

p = XI • 

k=l 



(57) 



(58) 



Closing the contour of integration in Eq. (57) on the left, we obtain in the notation 
of Ref. [38]: 



El^Q{m ,M ,a,P,ai,--- , cTg-^, pi, • • • ,Pr 



•1— Tl - 

I 7r2 



(-M2)2 



q—a—l3—(j—p 



r(a)r(^)r (§ ) 



jj r(f - CTj) 



n 



r(§-pfe)l r(c7-f(g-r-l)) 



4 r(^.) ( Ui^i r(Pfe) j r(f(g-r)-a) 

r(/3)r(§(g-r)-(7-;9)r(p + ^ + a-|(g-l))r(a + p + ;3 + a-fg) 

r(a + /3-f(g-r-l)) 

+ f (g- l),a + /3-Fa-Fp- §g 

M2 



(7 + /3-§(g-r -!),! + (7 + ;3-§(g-r) 



/ m2 \ /3 r (f (g-r)-a) T (cx+/3- f (g-r)) T (a+p-f r) T (p-f (r-1)) 

VM^; r(§) 



X 3^2 



/9-§(r-l),a+p-§r,f(g-r)-c7 
n,l + f(g-r)-c7-A 



M2 



(59) 
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Let us discuss the differential reduction of the two hypergeometric functions in 

Eq. (59) assuming that /3, aj (cj > 2), and pk are integers. We distinguish between 
two cases: r = 1 (g > 3) and r > 2 {q > 4). For r = 1, both hypergeometric 
functions in Eq. (59) are reducible to 2^1 functions with one unit upper parameter, 
namely 



1, 

/2-§(g-2) 



/2 + t 



(60) 



respectively. Standard approaches yield one three-loop bubble master integral and 
integrals expressible in terms of Gamma functions. For r > 2, the first hypergeomet- 
ric function is expressible in terms of a 3F2 function with one unit upper parameter, 
and the second one is reducible to a 2F1 function with both upper parameters con- 
taining non-zero e parts, namely 



(1,^)X3F2 



l,Ii-§(g-l),l2-fg 
/3-f(«-r-l),/4-§(g-r) 



Ii-§(r-l),l2- 

/3 + § 



(61) 

respectively. For q = 3 and r = 1, we present the explicit form of the master integral, 
with the powers of propagator being all equal to unity, in dimension n = A — 2e. It 
reads: 



2-£l 



(1-g) r(i+3g) 
3(i-3£) r(i+£) 



l,-l + 3£ 

1+s 



(M 
2 



2,i_3,r2(i-g)r(i + 2e)r(i + £) 



2e3(i-2£)(l-£) 



m 
M2 



+ 



^2^i-2e r(l+e) 



M2 



T{l-e) 



2F1 



l,e 
2-e 




(62) 



The two Gauss hypergeometric functions in Eq. (62) can be reduced to the basis 
functions considered in Refs. [16,20,35] by the following relations: 



(1 - 2e) 2F1 
(1 - 2e) 2F1 



1,8 

2-e 



l-e-£{l-z) 2F1 



1,1 + e 
2-e 



1,-1 + 3e 
1 + e 

3e{l - 3e)z{l 



l-2e-z{l- Se) 



+ 



1 + e 



2F1 



l,l + 3e 
2 + e 



(63) 
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For illustration, we present here the first few coefficients of the e expansion of 
Eq. (62): 

Ef20(^^M^ 1,1, 1,1,1) ^ .^2^1-3^ ( ^ + '^^ , 3 + z[5-31nz)] 

+ - I ?^±^ -{l-z) [Li2 (z) + ln(z) ln(l - z)] + z In z(ln z - 5) + C2 
e 

+(1 - z) [2Si,2(^) + 2 (Li2 (z) - C2) ln(l -z) + Lis (z) Inz] 
+(1 - z) {ln(l -z)lnz [Inz + ln(l - z)] - 5 [Lis (z) + ln(z) ln(l - z)]} 

45 + 49Z , ^ , 2 «/o , 5-2z^ 2 , 3 
H 172;ln2; + 52;ln^z + 2(3-zln2;)C2 5 — (3- ^zln"^ z 

+£{1 - z)i^U2 (z) [21n(l - z) (5 - ln(l -z)-\nz)- \nz{\nz- 5) - 4(2 - 17] 

+Li3 (z) Inz - Li4 (z) + 2Si,2(z) [5 - 21n(l - z) - Inz] - 4Si,3(z) 
+21n(l - z)C2 [ln(l - z) - 5 - Inz] + 41n(l - z)C3 



+ In z ln(l - z) [ln(l - z) + In z 
1 



2 

5 - -(ln(l - z) + lnz) 



+ lnzln(l — z 
1 



Inzln(l-z)-17 



+£<( zlnz { ^ In^ z(lnz - 10) + 171nz + 2C2(lnz - 5) + 2(3 - 49 

+^ + 43z + 25C2 - y (3 - ^)C3 + ^C4(23 - 14z) I + 0{e) ] , (64) 



where 



2 

m 



-SI- (M) 

To check Eq. (64), we evaluate the first few coefficients of the e expansion of the 
original diagram in the large-mass limit [42] using the program packages developed 
in Refs. [43,44] to find agreement. 

3.3 J22 

Let us consider the g-loop sunset diagram in Fig. 2 with two massive lines of 
the same mass m and q — 2 massless subloops, which is one of the most frequently 
studied Feynman diagrams. It is defined as 



[fcf j'^i • • • [A:^_i]'^'3-i [A;2 - [(ki + ^2 H \- kg + p)"^ - m2]«2 
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The MeUin-Barnes representation of Eq. (66) reads: 

(_^2)f-ai.2(^2)^(.-l)-a r(| -CTfc) 



[ii-"W2]-'?r(ai)r(a2) 



n 

.A;=l 



r(fTfc) 



X / (i^ 



p2 ^ * r(ai + i)r(a2 + t)r(ai,2 - f + t)r(§ + t)r(a - |(g - 1) - i) 



mV r(ai,2 + 2t)r(fg-c7 + i) 

where a is defined by Eq. (46) with s = q — 1 and 

"1,2 = ai + 02 • 



,(67) 



(68) 



Closing the contour of integration in Eq. (67) on the left, we obtain in the notation 
of Ref. [38]: 



,2 2 



22l,m ,p ,Q;i,Q;2,cri, 



o-<7-l 



r(ai)r(a2) \J-^i r(afe) 



r (ai + a - f (g - 1)) r (a2 + <7 - §(9 - 1)) r (^7 - f (g - 2)) T (ai,2 + a - f g) 



r(ai,2 + 2c7-n(g-l))r(f) 



Ctl + <T - f (g - 1), 02 + o- - |(g - 1) 



(7 



§, 5(ai,2 - n{q - 1)) + a, i(l + 0:1,2 - n{q 



2), 01,2 +<T - f ( 
-1)) + ^, 



Am? 



(69) 



For g = 1 (see Footnote 10), the 4F3 function in Eq. (69) is reduced to a 3F2 function, 
in agreement with Ref. [38]. In the two- loop case [q = 2), the hypergeometric 
representation of this diagram was derived in Refs. [45,46]. 

Let us analyze the reduction of the hypergeometric function in Eq. (69) assuming 
that all parameters, oi, 02; and ak, are integer. For q = 1 {ok = 0), we obtain a 
3F2 function with integer differences between upper and lower parameters, which, 
according to Criterion 1, is reducible to a 2F1 function with one integer upper 
parameter, so that wc have a basis hypergeometric function plus a rational function. 
(For the one-loop propagator, there are one master integral and bubble integrals.) 
For = 2, we get a 4F3 function with integer parameter differences and one integer 
upper parameter, so that it is reducible to a 3F2 function with one integer upper 
parameter and its first derivative. 



(1,^) X 3i^2 



n J 

2, -'2 



n 



(70) 



plus a rational function. In this case, there are two nontrivial master integrals of 
the same topology and bubble integrals, in accordance with the results of Ref. [47]. 
For g > 3, we have a 4F3 function with integer differences of parameters, which, 
according to Criterion 1, is reducible to a 3F2 function and its first two derivatives. 



{1,6,6' 



iF2 



/l-f(g-l),/2-i(g-2),/3-§^ 
§,/4 + ^-§(g-l) 



(71) 



IOt 



°In the one- loop case {q = 1), the factor |nfe=i 
Eq. (67) agrees with the one in Ref. [38]. 



r(<Tfc) 



■ I is equal to unity, and the representation of 
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3.4 -£^1220 ^1220 

Let us consider the g-loop bubble diagram -E^220 ^^S- ^ ^^^'^ massive lines of 
different masses and s plus r massless propagators, 

-^122o("^^^^«l'«2,/3,0■l,••• ,C7s,pi,--- ,Pr) = 

dJ'iki ■■■kg) 



I 



[kfjTl . . . [k2]as [A;2^^ _ M^]"! [{ki + ^2 + • • • + A:,+i + kg)^ - M^]^^ 

[k',+2]'' ■ ■ ■ [kl_^]Pr-r + ■■■ Vl + WY^ikl - ' ^^^^ 



where, by construction, s > and r > 2 and, as a consequence, g = s + r + l>3. 
The Mellin-Barnes representation of Eq. (72) reads: 

£;W-^M^«„«2,Aal,•••,-.,Pl,•••,P.) = r(a,)r(a2)r(^)r (|) 

r(f - 1 / fr r(f - pfc) 1 r (p - §(r - 1)) 



X dt 



\ ' r(ai + ^)r(a2 + t)r(ai,2 - f + t)r(f + t) 
My r(ai,2 + 2i) 

r (g - §s - r (g + /3 + p - f (g - 1) - t) r (§(g - 1) - a - p + t) 

r(f(s + l)-a + i) 



(73) 



where a, p, and ai_2 are defined in Eqs. (46), (58), and (68), respectively. 

Diagram -6^220 Fig. 2 emerges from £^^220 with r = 1 and s = g — 2 in the 

smooth limit p = pi — t- 0. Then, Eq. (73) simplifies due to p^ft^ ^j^^ — = 1 • 

Closing the contour of integration in Eq. (73) on the left, we obtain in the 
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notation of Ref . [38] : 



^1220 (m^, , ai , q;2 , /3, (Ti , ■ 



,Pi, 



rg— 01,2— o- 



r(ai)r(a2)r(/3)r(§) 




r(f-/,fc)lr(p-§(r-i)) 



fe=i 



5F4 



+ 



s'-^-" r (g - §(8 - 1)) r (f r - p) r (p + - §r) 
M2y) r (f ) r (ai,2 + 2a - ns) 

xr (ai,, - f r {a2,a - P) r (ai,2,. - §(s + 1)) 

ai,a - f S, a2,a - f S, Oi,2,a - f (s + 1), fJ - f(s - 1), §r - 
§ , + a, "^'^V""' + tr, f r - p + 1 - ^ 

r(/3)r(§r-aff,,)r(«.,^,,-^(g-2)) 
r (ai,2 + 2c7 + 2/3 + 2p - n(g - 1)) F [a^^p - §(r - 1)) 

xr {ai^a,l3,p - f(g - 1)) r (a2,a,/3,p - f (9 - 1)) r (ai,2,a,/3,p " 



4M2 y 







ai,o-,/3,p — §(<?-!), a2,o-,/3,p-§(g- 1), «i,2,(7,/3,p — f <?, ao-,/3,p — 2), /3 

1 + «/3,P -2^^ «<^,/3,P + 2 ' «<^>/3,P + 2 ' «/3,P " 2 " 1) 



m, 

4M2 



(74) 



where we have introduced the short-hand notations 

= ai + o- , a/3,p = ^ + p , 
ai,2,CT = ai + 02 + o- , a(T,/3,p = + ^ + P , 
ai,2,(T,/3 = ai + a2 + u + fi , ai,a,i3,p = Qi + a + f} + p , 
ai,2,a,i3,p = ai + a2 + a + P + p , i = l,2. 



(75) 



Let us analyze the result of the differential reduction of the two hypergeometric 

functions in Eq. (74) assuming that all parameters, ai, 02, /3, (Jk, Pj, sue integer. We 
distinguish between the three cases: s = 0, s = 1, and s > 2. For s = (a^ = 0, 
r = q — 1), both hypergeometric functions in Eq. (74) are reducible to 2-^1 functions 
with one unit upper parameter, 



2F1 



2F1 



(76) 



Standard approaches yield one master integral. For s = 1 {r = q — 2), each of the 
two hypergeometric functions is reducible to a 3F2 function with one unit upper 
parameter. 



(1,^)X3F2 



l,/i-n,/2-§ 



(1,^)X3F2 



l,/l-f((?-l),/2-§g 

/3-§(g-3),/4+i-§(g-i; 



(77) 

For s > 2, the first hypergeometric function is reducible to a 3F2 function and its 
first two derivatives, while the second one is expressible in terms of a 4F3 function 
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and its first two derivatives, namely 



{1,9,9) X,F,[ l^^n j^+l_n, 

n««2. l,/i-fg,/2-f(g-l),/3-f(g-2) 



(78) 



respectively. Notice that, according to Criterion 4, a third derivative does not 
appear because one of the upper parameters is integer. 

We now present an explicit result for the three-loop case in which there is no 
massless propagator inside the massive loop, q = 3, s = 0, = a = 0, r = 2, and 
n = 4— 2e. In this case, there is one master integral with ai = q;2 = /3 = Pi = P2 = 1 
("1,2 = P = 2), 

l3 / ,^^xi_3£ 



E?22o(m2, M2, 1, 1, 1, 1, 1) = [iTT^-^] (M2)1- 

Tjl - e)r\l + e)Til + 2e) i fm^y-'' (l,e 



m 



4M2 J 



r(l + 2e)r(l + 3e) 



3(1 -3£) r(i + 4£)r(i + £) 



1,-1 +3e 



m 

4M2 



(79) 



To reduce the Gauss hypergeometric functions in Eq. (79) to sets of functions studied 
in Refs. [16,20,35], we apply the following relations: 



(l-2£) 2^1 (V 

/l,-l + 3e 

U2£ 



z]=l-2e{l-z) 2Fi(^'\+' 

^ , (l-3g) , 12£(1 - 3g)z(l - z) /l,l + 3e 
-■^ ^^l_2e+ (i_2e)(l + 4e) ''^^l^| + 2e 



(80) 



The first few coefficients of the e expansion of Eq. (79) are given by 



£;f22o(m^M^ 1,1, 1,1,1) 
r3(l + e) [iir^-'f 



(M 



2\l-3e I 



(1 + z) 6 + z(5-31nz) 



3£3 



+ 



3£2 



+— {25 + z(17- 151nz + 31n2z + 3C2)} 



+2 



(4-^) ji-y) 
e (1 + 2/) 



Li2 (1 - y) + ^ In^ y + ^ In z In y 



+^(90 + 49z) + ^C3(8-7z)-z 



171n2;- 51n^2; + -ln^2; 



+ C22(5-21nz) 



1 + y 



{A-z)\ 4Si,2(y) + 2Si,2(y') - 4Si,2(-y) - GLia (y) - 2Li3 (-y) 



+Li2 (-y) [4 ln(l - y) - 2 In z] + Li2 (1 - y) [10 - 8 ln(l - y) - 4 ln(l + y)] 
+ ln(l -y) [2C2 -41nyln(l -y) -31n2y] +ln2y 



5 2 
--ln(l+y) + -lny 



-C2(ln^ - 41ny) - (3 + ln2;lny 



5 + - Iny - 21n(l + y) - In 2; 



+ 0{e)\ ,(81) 
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where z is defined in Eq. (65) and the variable y is defined as 



z-4 



1 + 



(1 - y? 
y 



4-z 



y 



(82) 



The 0{e) term is too long to be presented here, but is available from Ref. [48]. We 
mention here only that, in accordance with the expansion of the Gauss hypergeo- 
metric function constructed in Refs. [20,35], this term is expressible just in terms of 
Nielsen polylogarithms. To check Eq. (81), we evaluate the first few coefficients of 
the £ expansion of the original diagram in the large-mass limit [42] using the pro- 
gram packages developed in Refs. [43,44] to find agreement. In order to construct 
the series expansion about 1 — y, which serves as the small parameter of the large- 
mass expansion, we employ the trick described in Ref. [49]. In fact, the variable 
1 — y can be written as a series in the small parameter z as 



tVz 



1 



z z^ {2ky. 1 



8^^^(fc!)2fc + 



1 V16/ 



■ (83) 



For completeness, we also present the explicit result for diagram S^22o i^i Fig. 2:^^ 



(-M2)S«-"i.2-cT-/3p^„_^^ 



■q-2 

n 



r(§-afc) 



r(ai)r(a2)r(/3)r(f) 

/ r (ai,,,^ - |(g - 1)) r (a2,.,;3 - f (^ - 1)) T (ai,2A. - f g) r (a,,^ - f (g - 2)) 



r (ai,2 + 2/3 + 2c7 - n(g - 1)) 



«1,/3,C7 - §(g - l),02,/3,o- 
l(«l,2 



/ vn? 



§(g-l),ai,2 + /3 + a-fg,a<,,^-§(g-2] 

n{q - 1)) + a + /3, i(ai,2 + 1 - n(g - 1)) + a + /3, 1 + /3 - § 

r {ai,a - f (g - 2)) r (q2,. - f (g - 2)) r (ai,2,. - f (g - 1)) 
... - ,(^_2)) 



X 



r (f - /3) r (ai,2 + 2a 
r(a-f(g-3))r(/3-§) 

(g - 2), 02,. - f (g - 2), ai,2,a - f (g - 1), a - f (g - 3) 




/ ai,a - 2 W - a2,<7 - 2 w - ^J, ai,2,<7 - 2 W - J-j, o" - 2 W ~ 
^ ^ V ^(«i,2 - "(g - 2)) + CT, i(ai,2 + 1 - n(g - 2)) + a, 1 - ;3 + § 

where g > 2. For g = 2, the hyper geometric representation was derived in Ref. [51].^^ 
The result of the differential reduction of the hyper geometric functions in Eq. (84) 
assuming that all parameters, ai, a2, /3, and cr^, are integer may be derived for 
g = 2, g = 3, and g > 4 by directly substituting r = 1 and s = g — 2 in Eqs. (76)— 
(78), respectively. We only point out here that, for g = 3, the second hypergeometric 



^-'^For the four-loop case {q ~ 4) with equal masses m — M, the Mellin-Barnes representation of this 
diagram was published in Ref. [50]. In order to recover the result of Ref. [50] from Eq. (73), it is sufficient 



to redefine the variable of integration t tohe t = N/2 — 
^^In this case, we have a = and 1 11^=1 ~T^o^^} ~ 



Q!12 

1. 



s. 
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function in Eq. (77) is reducible to a 2F1 function and its first derivative, so that 
the differential reduction yields 



(l,^)x 2F1 



/i-n,/2-f 
^ + /3 - n 



(1,^) X 3F2 



l,Ji-n,/2-i 



z] , (85) 



and that, for g > 4, the second hypergeometric function in Eq. (78) is reducible to 
a 3F2 function and its first two derivatives, so that the differential reduction yields 



(1,9,9^) XsF2 



Ii-fg,/2-t(g-l),/3-f(g-: 
i+/4-f(g-i),/5-i 

h - §(g - 2), 72 - §(g - i),/3 - f (g - 3) 
/4 + f,/5 + ^-f(g-2) 



(86) 



It is interesting to note that, in the single-scale case m = M [52], there is only 
one master integral [43]. This is because the hypergeometric functions 3^2 and 
2-Fi with special values of parameters and argument z = 1/4 are expressible as 
products of Gamma functions. For details, see Eqs. (4.36) and (4.42) in Ref. [29]. A 
diagrammatic interpretation of similar identities was presented in Ref. [53] . On the 
other hand, this is a manifestation of the existence of functional relations between 
Feynman diagrams [54]. 

We now present the explicit result for the three-loop case g = 3, n = 4 — 2e. In 
this case, the first master integral corresponds to ai = a2 = /3 = cri = 1, while we 
choose ai = P = a\ = 1 and a2 = 2 for the second one. 



B 



1220 



(m2,M2,l,l,l,l) = (M^ 



2\2-3£ 



r3(l + £)(m^-^)3 
e3(l-£)2(i_2£) 



2(1 - £)(1 - 4£)r^ (1 + 2e) r (1 + 3£) r(l - £) 
3(1 - 3£)(2 - 3£)r(l + 4£)r2(l + e) 



2F1 



-1 + 2£,-2-F3£ 
-i + 2£ 



4M2 J 



+ 



l-e 



iF2 



1,£, -1 +2£ 

2-e,^ + e 



m 



4M2 J 



Bl^^^{m\M\l,2,l,l) = (M 



2£3(1 



2(1 - e)(l - 4£)r2 (1 + 2£) r (1 + 3£) r(i - £) 
3(1 - 2e)(l - 3£)r(l + 4£)r2(l + e) 



iFi 



-l-h2£, -l + 3£ 

-i + 2£ 



m 

4M2 



+ 



WJ 



l-e 



iF2 



1,£,2£ 
2-£, i-F£ 



m 



4M2 J 



(87) 



In order to express all hypergeometric functions entering Eq. (87) in terms of func- 
tions which were analyzed in Ref. [20] (see Appendix A for details), we apply the 
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following set of relations: 



24:z{l- z){l + 2z)s^ /l + 2e,l + 3e 



(1 + 4£) 



+ [1 - 4£ - 4z + 8ez^ + Uez] 2^1 2^ ^) 



24z(l - z)e^ „ n + 2e, 1 + 3e 

2-^*1 



1 +4£ 



I + 26 



+ [1 - 4£ - 2z + lOze)] 2F1 (I 

V 2 

2z(l-2£)(l -3e) 



2e,3£ 



l-£ 



^^^t 2-£,i+£ 



l-2£ 



+ [2£(1 - Az) - (1 - 2z)] 3F2 ^ ^ 
.(l-3£)(2-3£)3F2(^2-7,'iV£ 



8z(l - ^)£2 
+ - 2-£,i+£ 



(l-£)(l + 2£) 



-(l-£)(l-2£ + 2z£) 



2, 1 + e, 1 + 2£ 
3 
2 



A£^z{l- z){l + 2z) 



3-^2 



2,l + £,l + 2e 



(1 + 2£) ^ V 2 - £, § + e 

[l-Az- 2£(1 - 8. - 2^2)] 3i^2 ( 1 



1,£,2£ 

+ £ 



' 2 



Here, the following relation was used: 

z{h - 02) (6 -az){l + h- ai){f - 03) 3F2 
= -6(1 - az){l -z)[a^-f + (02 - h)z] B 3F2 
+6(1 - aA a2{a2 - h)z^ + (03 - /)(1 - /) - 6(1 + 02 + 2a3 - h)z 



1,02,03 - 1 
1,02,03 



- [/(2a2 + 03 - 26) + a2(o2 - 1) - (02 + 03)^] z| 3i^2 ( ^' ^) 



-6(/ - 1) [(1 - a3)(/ - 03) - (6 - 1)(6 - a2)z] . 



(89) 
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The first coefficients of the £ expansions of Eq. (87) read: 



Bl 



1220 V"*- 5 



.M2, 1,1, 1,1) 



r3(l + e)(i7r2-^)3 

L r 25 20 
e I 12^y 



(M^) 



2-3e / 1 + 2z ^ 1 f 7 



3e3 ' £2 ]6 ' 3^ 12'' 



^ 2; In 2; 



1 r 25 20 5 2 1 , r 



8.,-, ^ 5 35 145 2 1,3 (16-z) , 2 (80 - 15z) , 



-4(1 - z) [Si,2(l -y) + InyLis (1 - y)] - {I - z) In^ 



^l'220K.M',l>2,l,l) 

r3(l + e)(f7r2-e)3 



Li2 (1 - y) + ^ ln2;Iny + ^ In^ ; 



In z + - In y 
o 



+ 0{e)\ , 



= (^^)"'M^ + ^{4 + 5.-3.In4 



+7^ {1+42;} - -^zln2;{4- Inz} - ^ + ^z-^zlnz 
6e 4e 3 b 2 



1 — In 2 + - In z 



+(4 - z, 



(1-^) 
(1 + ?/) 



2Li2 (1 - y) + ^ In^ y + In 2 In y 



+2(z-2) 



Si,2 (1 - y) + In yLi2 (1 - y) + ^ In^ y + i In z In^ y - ^Cs 



+ 0{e) , (90) 



where y and 2 are defined in Eqs. (82) and (65), respectively. The 0{s) term is too 
long to be presented here, but is available from Ref. [48]. To check Eq. (90), we 
evaluate the first few coefficients of the e expansion of the original diagram in the 
large-mass limit [42] using the program packages developed in Refs. [43,44] to find 
agreement. 



3.5 Vi^220 and 

•^1220 

Let us consider the g-loop on-shell propagator diagram F]^220 ^^S- ^ ^^^'^ * P^^^ ^ 
massless propagators, 

^122o("^^'*^^''^l'«2,/3,0-l,--- ,as,Pl,--- ,Pr) 

(P{kl ■■■kg) 



[fc2]<Ti . . . [^2]a. [/,2^^ _ M2]"l [{kl +k2 + --- + ks+1 + kqf - M2]"2 



[kl^2Y' ■ ■ ■ [kl-iY-^ [{ks+2 + ■■■ + k,-i + k,)-')]Pr {{k, - V? - 



,(91) 



27 



where q = s + r + 1. The MeUin-Barnes representation of Eq. (91) reads: 

(-M2)t~°i'2(_77j,2)§(9-l)-o— 
^122o('^^-^^"l'<^2,^,Crl,••• ,Crs,pi,--- ,pr) = 



1 "1 



r(Qi)r(a2)r(/3) 

n-i g 



n ^^t-^.) jfr rd-Pfc) r(p-§(r-i)) 
f=i (iM / r(tr-,) 



X \ dt\—r\ 



m2 \ * r(ai + t)V{a2 + t)r(ai,2 - f + t)V{\ + i) 



r(ai,2 + 2t) 

r (ct - §s - i) r (a + /3 + p - §(g - 1) - r (n(^ - l) - 2(j - ^ - 2/9 + 2t) 
r (§(s + 1) - a + t) r (^g - a - /? - /9 + t) 



(92) 



where a, p, and qi^2 are defined in Eqs. (46), (58), and (68), respectively. 

Diagram J^na -^iS- ^ emerges from ^^220 ^^^^ r = 1 and s = g— 2 in the smooth 
hmit /9 = pi ^ 0. For the special case g = 3, the Mellin-Barnes representation has 
recently been presented in Rcf. [55]. 

Closing the contour of integration in Eq. (92) on the left, we obtain: 



n \ q 

^l22o("^^^^«l'«2,/3,o■l,••• ,0-5, pi,-- - ,p 



( n \q ^ n 



r(ai)r(a2)r(/3)r (§ ) 




j^r(§-p,)lr(p-f(r-i)) 



r(Pfc) j r(§r-p) 



X < ^ 



m2 \ t"-"''''' r (a - § (s - 1)) r (nr - 2p - /3) r (p + /3 - f r) 



M^J r(ai,2 + 2a-ns)r(f(r + l)-p-/3) 

X r (ai,, - §s) r (a2,, - f s) r (ai,2,. - f (s + 1)) 



2'5,a2,a- 2'5)ai,2,<7- 2('5+l)><7-2('S-l), ^^-P, - P 

y §, fiif^ + cx, |(r+l)-/3-p, f r-^-p+1 

r(fr-a^,,)r(a.,^,p-f(g-2))r(/3) 
r (ai,2 + 2c7 + 2/3 + 2p - n{q - 1)) T (0^3,^ - f (r - 1)) 



X 



+ 



X r (ai,^,/3,p - f (g - 1)) r (a2,a,/3,p - f (g - 1)) r (ai,2,a,/3 
X 6F5 



ai,o-,/3,p-§(g'-l),a2,<T,/3,p-§(9-l),ai_2,<T,/3,p-f9,«(T,/3,p-§(^-2) - — 



_ ' 2 ' 2 

n 1, ^ n„ ai,2-»('?-l) , „ Qi,2+l-»(<?-l) , „ n / -, n 

2) -L+a/3,p-2^ 2 I-"ct,/3,p, 2 \-a<T,/3,p,ai3,p-2{^-^) 



m 
M2 



where ai,o-, a2,a, a^^p, ai,2,<7, aa,/3,p, «i,<t,^,p) and a2,a,/3,p are defined in Eq. (75). In the 
special case in which there is no massless propagator inside the massive loop, s = 

(r = g — 1, cr = 0, 11^=1 rVfc) ~ 6-^5 functions in Eq. (93) are reduced to 

5F4 functions. For g = 2, s = 0, and r = 1, the hypergeometric representation was 
given in Ref. [56]. 

Let us analyze the result of the differential reduction of the hypergeometric 
functions in Eq. (93) assuming that all parameters, ai, a2, P, (Tk, and pk, are 
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integer. For r 7^ 0, we distinguish between the three cases s = 0, s = 1, and s > 2. 

For s = 0{r = q— l,a = 0), each of the two hypcrgcomctric functions in Eq. (93) is 
reducible to a 3F2 function with one unit upper parameter and its first derivative, 



i3 + 2! -'4 + -2? 
l,/l + i,/2-§g 



1: 



(1,^) X 3F2 



-^3 + §) 2 + -^4 



(94) 



For s = 1 (r = g — 2), each of the two hypergeometric functions in Eq. (93) is 
reducible to a 4F3 function with one unit upper parameter and its first two deriva- 
tives, 



{i,ey)x 4F3 



l,Ii-n,l2-f,/3 + i + §(g-2) 
l,/i + i,l2-§(?,l3 -§((?-!) 

/4 + §,i+/5-^(g-l),/6-§(?-3) 



(95) 



For s > 2, one of the hypergeometric functions in Eq. (93) is reducible to a 4-F3 
function with all upper parameters having non-zero e parts, whereas the second one 
is expressible in terms of a 5F4 function with one unit upper parameter, 



h - f (s - 1),/2 - %s,h - f (s + 1),/4 + \ 

/5 + f,/6 + ^-fs,/7 + f(r + l) 



h + ^,^+h-^{q-l),h 



(96) 



For completeness, we also present an explicit result for diagram J1220 ™ ^i§- 

•^122o("^^'-^^'"l-"2,,^,0-l, • • • ,0-9-2) 



ii-"7r2 



T(/3- §) (-M2)t9- 



■C«l,2,( 




r(«i)r(a2)r(^)r(f) T{ak) 

.2 X 2-/3 r (a-f (g-3)) r (ai,.-§(g-2)) F (02,.- f (g-2)) F (ai,2,.-f (g-1)) 



r(ai,2 + 2(7-n(g-2)) 



ai,.-f(g-2),a2,.-f(.Z-2),ai,2,.-§(«-l),a-t(<z-3),^,^ 



n 

2 ' ' 



n-/3,f + l-^, 
F (§-/?) F(a.,^-f(g-2))F(ai 



o:i,2--»(q^2) Q:i,yH-n(q^2) 



+ 



X 6i^5 



n 

<^,/3- 2 



2 ' 2 

((7-1)) F (a2,,,;3-§(g-l)) r K2,a,/3-i<?) 



F {p- f ) F (ai,2+2a+2^-n(g-l)) 



ai,a,/3-f (g-l),a2,^,/3-f (g-l),ai,2,^,/3-ig,a<T,/3-t(g-2) ^ ^ 



2 ' 2 



m 
M2 



(97) 



In this case, the result of the differential reduction assuming that all parameters, 
ai, a2, P, and cr^, are integer can be derived for g = 2, g = 3, and g > 4 by directly 
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substituting r = 1 and s = g — 2 in Eqs. (94)-(96), respectively. We only point out 

here that, for q = 3, the second hyper geometric function in Eq. (95) is reducible to 
a 3F2 function and its first two derivatives, so that the differential reduction yields 



{i,e,e^) X 4i^3 



l,/i-n,l2-f,/3 + i + § 
/4 + f ) 5 + ^ — f ) -^6 + ^ 



{1,9,9'^) X 3F2 \ r 7 n 1 



3n 

I4 + ^A + h-n 



n 1 
2 ' 2 



(98) 



and that, for g > 4, the 5F4 hypergeometric functions in Eq. (96) is reducible to a 
4i^3 function with all upper parameters having non-zero e parts: 



{1,9,9^6^) X 4F3 



A - f (g - 3), /2 - f (g - 2), /3 - f (g - 1), I4 + i + 



h + %,h + k-^{q-2),Ij + 



n 



{1,9,9^9^) X 4i^3 ' + ^'^2 - ^q,h - i(g - l),/4 - t(Q 



/5 + f,i+/6-f((?-l),/7 



(99) 



In particular, for g = 2, in accordance with Ref. [57], there are two nontrivial master 
integrals plus bubble integrals. For g = 3, in accordance with Ref. [55], there are 
three nontrivial master integrals plus bubble integrals. 



3.6 Two- loop vertex 

Let us consider the diagram F2 shown in Fig. 2 

F2{M^,p'^,ai,a2,ai,a2,) = 



I 



(P{kik2) 



[{ki - pifY^[{ki - p2fY^[klf[kl - M2]°i[(A;i - ^2)^ - M2]«2 

2 (-M2)i^i2(-p2) 1^12-/3 



pI=pI=q 



X / dt 



T{ai)T{a2)T{ai)T{a2) 

p2 y r(-t)r(ai+t)r(a2+t)r(ai2-f +t) 

M'^J r(ai2 + 2t) 

r (|-c7i-^+t) r {r^-a2-P+t) V {ai2+P-^-t) 
T{n -ai2-P + t) 



(100) 
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where a is defined in Eq. (46) with s = 2 and 0:12 in Eq. (68). Closing the contour 
of integration on the left we obtain in the notation of Ref. [38] : 



F2 (M^ , , ai , 02 , CTl , (72 , ) 



2 (-M2)2^12 

r(ai)r(a2) 

/r(ai)r(a2)r (f-ai-^) r (f -a2-^) r ((Ti2+^-f ) r («i2-f ) 



r(ai)r((72)r(ai2)r {n-au-P) r (f -(712-/3) 



+ 



AM^ J 

^\^'^2 r (ai+c7i2+;3-§) r (a2+ai2+;3-§) r(ai2+ai2+^ 



. Q;i,Q;2,ai2-f, §-0-1-/3, §-0-2-^ 

5^ ft ai2 ai2+l 

n — (712- P, 1+2 —(712 — P, — 



-n 



M2 J 



r(ai2+2ai2+2^-n)r(f) 



'c7i,c72,ai+<Ti2+/3-§,a2+cri2+/3-f ,ai2+(7i2+/3-n 



4M2 y 



(.101) 



By differential reduction, the two hypergeometric functions in Eq. (101) may be 
reduced to 



(i,^)x 3i^2 ' ! 



P 



4M2 



(1,^) X 3F2 



1, 1,/i-n 

-^2 + 1, ^3+5-5 



P 



4M2 



respectively, where {/„} is a set of integers. In accordance with Ref. [20], the first few 

coefficients of the e expansion of this diagram are expressible in terms of Remiddi- 
Vermaseren functions, i.e. multiple polylogarithms of the square roots of unity, of 
argument 

i-'-^M- (103) 



1 + 



z 

2+4 



where z = p'^/M'^ (sec also the discussion in Ref. [58]). The on-shell value z = 1 cor- 
responds to Remiddi-Vermaseren functions of argument (3 — v^)/2 = [(a/5 — l)/2]^, 
in agreement with results of Ref. [59] . It is interesting to note that these constants 
are generated not only in Higgs-boson decay [59] , but also in orthopositronium decay 
[60]. 



(102) 



4 Master integrals in differential reduction and 
IBP technique 

In this section, we compare the differential-reduction and IPB techniques with re- 
spect to the numbers of master integrals which they produce. For this purpose, 
we return to the Feynman diagrams depicted in Figs. 1 and 2. They all have the 
following hypergeometric structure: 

k 

$(n,j;z) = Y,^^''Clain,j)p+iFpi^a;Ba;Kz) , (104) 

a=l 
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Table 1: Highest powers of the differential operator 9 generated by the differential reduc- 
tion of the Feynman diagrams shown in Fig. 2. 



h 


77K? 

-^120 


s 


-^1220 


*^1220 


q 


•^22 


-"1220 


"^1220 


1 


1 





1 


2 


1 


1 






> 2 


2 


1 


2 


3 


2 


2 


1 


2 






> 2 


3 


4 


3 


3 


2 


3 












> 4 


3 


3 


4 



where j is the set of the powers of the propagators of the Feynman diagram, n is 
dimension of space-time, k is the number of hypergeometric functions, 
are linear combinations of j and n with rational coefficients, k is a rational number 
(being k = 1, 1/4 in the considered cases), and Ci^ are products of F functions with 
arguments only depending on n and j. Being a sum of holonomic functions, $(j; z) 
is also holonomic. Thus, the number of basis elements on the r.h.s. of Eq. (104) is 
equal to the number of master-integrals {^^{z)} that may be derived from the l.h.s. 
of Eq. (104) by applying the integration-by-parts technique, which may be written 
symbolically as 

h 

^(nj-z) = ^Bkin,j;z)^k{n;z) . (105) 
k=l 

Here, it is understood that diagrams that are expressible in terms of Gamma func- 
tions [61] are not counted. The number of basis elements in the framework of 
differential reduction is defined to be the highest power of the differential operator 
e in 

V 

p+^Fp{A;B;z) = ^Pi{z)9's+iFs{A-h;B-l2;z) , (106) 

where /i,/2 are sets of integers, Pi{z) are rational functions, and the differential- 
reduction algorithm for exceptional values of parameters of hypergeometric function, 
constructed in Sections (2.3) and (2.4) is employed, in which case one has v < p in 
general. Our analysis demonstrates that there is a very simple relation between the 
number h of nontrivial master integrals follows from IBP, which are not expressible 
in terms of Gamma functions, and the maximal power v oi 9 generated by the 
differential reduction in Eq. (107), namely 

h = v + l. (107) 

This relation does not depend on the number k of hypergeometric functions entering 
Eq. (104). In Table 1, these highest powers are collected for the Feynman diagrams 
shown in Fig. 2. 
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5 Discussion and conclusion 



The idea that any Feynman diagram can be associated with a generahzed Horn- 
type hypergeometric function or a hnear combination of such functions was born at 
the end of the 60s [7] and was apphed during the last twenty years [38,46,62]. The 
MelHn-Barnes representation (1) of a Feynman diagram is universal. Under certain 
conditions, it can be converted to a linear combination of Horn-type hypergeometric 
functions (2). The latter possess very interesting properties. In fact, the systems of 
differential equations which they satisfy are sufficient for (i) the differential reduction 
of the original functions to restricted sets of basis functions, whose number follows 
directly from the system of differential equations; (ii) the construction of the all- 
order £ expansions of the basis hypergeometric functions in form of iterated integrals 
[16-19]. To our knowledge, the first property has never been discussed in the context 
of a reduction procedure of Feynman diagrams. 

The aim of the present paper was to show how the differential-reduction al- 
gorithm can be applied to reduce Feynman diagrams without having to exploit 
the integration-by-parts technique [15] or its dimensional generalization [37]. In 
our analysis, we considered a phcnomcnologically interesting class of Feynman di- 
agrams, namely those shown in Figs. 1 and 2, which are expressible in terms of 
generalized hypergeometric functions p+iFp of one variable and cannot be treated 
with currently available program packages. 

As first steps in our analysis, Takayama's differential-reduction algorithm [22] 
was extended to the case of exceptional values of parameters in Section 2.3, and 
the basis of differential reduction was written explicitly in Section 2.4. In partic- 
ular, this algorithm allowed us to complete the proofs of the theorems regarding 
the construction of all-order e expansions of hypergeometric functions presented in 
Refs. [16-19]. For non-exceptional values of parameters, this algorithm is imple- 
mented in the Mathematica based program HYPERDIRE to be described in a in a 
separate communication [63]. 

The differential-reduction formalism via the construction of step-up and step- 
down operators can be applied to the reduction of any Feynman diagram. Its 
advantages reside in the simplicity and universality of the construction of this type 
of operators, the full control over the analytical structure of the Feynman diagrams 
via the hypergeometric representations, the simplicity of the criterion of rcducibility 
of hypergeometric functions to simpler functions, and the existence of a few, recently 
developed algorithms [16-20,27] for analytically constructing the coefficients of the 
£ expansions of hypergeometric functions. 

The main criteria of rcducibility of hypergeometric functions, namely integer 
difference between upper and lower parameters and integer value of one of the 
upper parameters, are also valid for Horn-type hypergeometric functions of more 
than one variables. Moreover, the first type of reduction can be applied directly to 
Mellin-Barnes representations of Feynman diagrams. 

This approach becomes incomplete when some of the variables of a Horn-type 
hypergeometric function belong to its surface of singularities; for hypergeometric 
functions of one variable, this corresponds to z = 1. In this case, the extended Ore 
algebra technique [64] , the telescoping approach [65] , or the Laporta algorithm [66] 
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can be employed to obtain exact analytical results. In addition to these algorithms, 
there is another approach based on the analytical structure of the coefficients of the 
£ expansions of hypergeometric functions: instead of looking for exact recurrence 
relations, the limit z — > 1 can be taken for the coefficients of the Laurent expansions. 
In particular, this procedure can be applied to any hypergeometric function whose e 
expansion has coefficients that are expressible in terms of multiple polylogarithms. 
This algorithm may be applied to evaluate Feynman diagrams, if just the first few 
coefficients of their e expansions are required. 

As an illustration of the usefulness of our differential-reduction procedure, we 
considered a few examples, namely the Feynman diagrams shown in Figs. 1 and 
2 with arbitrary powers of propagators and space-time dimension, and explained 
how to express them in terms of hypergeometric functions and to construct their e 
expansions to higher orders. 

An interesting observation we made here is that, up to products of one- loop 
bubbles, the number of master integrals constructed through the integration-by- 
parts technique is equal to the highest power of the differential operator 6 generated 
by the differential reduction of hypergeometric functions plus one, independently of 
the number of hypergeometric functions occuring in the expressions for the Feynman 
diagrams. 
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A Sets of basis functions for the Laurent ex- 
pansions of hypergeometric functions 

A.l All parameters are integer 

Using the algorithm described in Ref. [18], the derivatives of a basis hypergeomet- 
ric function with integer parameters may be expressed in terms of the p^pj^^\z) 
functions defined in that reference, as^^ 

(p \ oo 
n ^''pp+k~'\^) ' i = 0, • • • ,p - 1 . 
r=l / fc=0 

(108) 

In Eq. (2.11a) of Ref. [18], the last parameter should have the value p — 2 rather than p —1. 
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The iterative solution of the p^^^~^\z) functions is given by Eqs. (2.13a)-(2.13b) 
in Rcf. [18], and their cxphcit forms in terms of generahzed polylogarithms are 
cohected in Section 3 therein. 



A. 2 One lower parameter is half- integer 

In Refs. [20,29], the e expansions of the functions 



p+iFp I 



{l + a,£}^,{2 + d,£}^-^+i 
V2 + /£,{! + ej£}^, {2 + Q£} 



P-R-l 



(109) 



with K ~ R > 2, were constructed up to functions of weight 4 and, in the case of 
z = 1/4, up to weight 5.^"^ For these parameters, after some trivial factorization, the 
coefficients of the e expansions contain functions of only one definite weight. In this 
section, we present explicit relations between the basis functions of the differential 
reduction, defined by Eq. (17), and the functions of Eq. (109). One of the auxiliary 
relations intensively used in the following reads: 



. 01,02, • • • ,ap 



l + z 



1, 1 + ai, 1 + 02, ■ • • , 1 + Op 



A.2.1 2F1 

A detailed discussion of the 2-^1 function was presented in Refs. [16,19,35]. 
A. 2. 2 3F2 

For the hypergeometric function 

' h + aie, I2 + a2e, h + a^e 



3-f^2 



(111) 



where {7^} are arbitrary integers and all Oj 7^ 0, the basis of the differential reduction 



IS 



oie, 026, 036 



2 + /£, 1 + Ci£ 

For P = 2, there are four functions of the type of Eq. (109): 



(112) 



1 + oie, 1 + 02£, 1 + ose, 1 

| + /£,2 + Cl£,2 



, 1 + ai£, 1 + 02£, 2 + 03£ 
' 2' | + /£,2 + Cl£ 



.1 + Oi£, 1 + 02£, 1 + 03£ 
|+/£,2 + Cl£ 

,1 + Oi£, 1 + 02£, 1 + 03£ 
^ |+/£,1 + Cl£ 



Z .(113) 



14 



The series representations of these hypergeometric functions are proportional to J^jLo jH-n-i 
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The functions of Eqs. (112) and (113) are related as 



3F2 



aiE, a2£, as£ 
=^^2 I 1 + fe, 1 + cie 



ai£, a2£, ciss 
J + /£, 1 + Cl£ 



1 + 



2z 



2aia2<i:i 



-3 



(l + 2/£)(l + Cl£) 

.3 



(l + 2/£)(l + Ci£) 



4-P3 



5-^2 



1 + OiS, 1 + 026, 1 + 036, 1 
| + /£,2 + Ci£,2 



1 + Oi£, 1 + a2£, 1 + 036 
I + /£, 2 + CIS 



(114) 

The second derivative in Eq. (112) may be obtained from any of the following 
relations: 



3-^2 
3i^2 



where 



1 + oi, 1 + 02, 1 + 03 
l + /,ci 

1 + oi, 1 + 02, 2 + 03 
1 + /, 1 + ci 

o,: Oi£, 



l/g(g + ci-l) 
z 010203 
1 cif9{e + a3) 
2010203(1 + 03) 



3-^2 

3i^2 



01,02,03 

/, Cl 
01,02,03 
/, Cl 



Cl — )■ 1 + Cl£, 



(115) 



(116) 



are to be substituted. In particular, it can be expressed as a linear combination of 
the second, third, and fourth functions in Eq. (113) as 



(71+72) fci 2 p 



01,02,03 

Cl,/ 



= 71 Cl 3-^2 
- [71 (ci - 1) + 7203] 3-P2 



Z +72(1 + 03) 3-^2 



1 + oi, 1 + 02,2 + 03 
l + ci,l + / 



(117) 



01O2O3 

1 + Ol, 1 + 02, 1 + 03 
ci,l + / 

1 + Ol, 1 + 02, 1 + 03 

1 + Cl, 1 + / 

where 71 and 72 are arbitrary numbers. Putting 72 = —71, we obtain a linear 
relation between the hypergeometric functions on the r.h.s. of Eq. (117). We checked 
that the £ expansions of the hypergeometric functions constructed in Ref. [20] and 
collected in Ref. [48] satisfy this relation. 

In the case oi = 1, the basis of the differential reduction is 

1,1 + 02£, 1 + 03£ 
| + /£,2 + Ci£ 

The first derivative in Eq. (118) can be obtained from any of the two functions 



{1,0} X 3F2 



(118) 



1, 1 + 026, 1 + 036 
3 
2 



3 + /£, 1 + cie 



iF2 



1, 1 + 02£, 2 + ase 
3 
2 



3 +/£,2 + Cl£ 



(119) 



by the single application of the step- up and step-down operators defined in Eq. (7). 
The result is 

1,1 + 02£, 1 + 03£ 



I + /£, 2 + Cie 



(71 + 12)0 3F2 

71 (1 + 03£) 3F2 
- [71(1 + 03£) + 72(1 + Cl£)] 3F2 



1, 1 + 026, 2 + 036 
3 
2 



I + /e,2 + cie 



z + 72(1 + ci£) 3-^2 



1,1 + 026, 1 + 036 
I + /£, 1 + Cie 



1,1 + 02e, 1 + 03e 

§ + /£, 2 + Ci£ 



(120) 
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where 71 and 72 are arbitrary numbers. 



A.2.3 4F3 

Similarly, for the hypergeometric function 

' Ii + aiE, I2 + a2£, I3 + a^e, I4 + a^e 



4F3 



(121) 



where {Ik} are arbitrary integers and all / 0, the basis of the differential reduction 
is 



{i,e,e\e^}x 4F3 ( 1 



ai£, a^E, asE, a^E 



5 + /£, 1 + CiE, 1 + C2E 



(122) 



There are seven functions of the type of Eq. (109), namely 



5-^4 
4-^3 
4i^3 
4i^3 



1 + aie, 1 + 026, 1 + a^e, 1 + a^e, 1 
I + /e,2 + cie,2 + C2e,2 

1 + ai£, 1 + 026, 1 + 036, 2 + 046 
I + /£, 2 + Ci£, 2 + C2£ 

1 + ai£, 1 + a2£, 2 + a^E, 2 + a4£ 

I + /£, 2 + Ci£, 2 + C2£ 

1 + oi£, 1 + a2£, 1 + a^e, 1 + a^e 



4-F3 

4i^3 
4i^3 



1 + ai£, 1 + 026, 1 + a3£, 1 + a^e 

I + /£, 2 + Ci£, 2 + C2£ 

1 + ai£, 1 + a2£, 1 + a3£, 1 + a4£ 

I + /£, 2 + Ci£, 1 + C2£ 

1 + ai£, 1 + a2£, 1 + a3£, 2 + a4£ 

I + /£, 1 + Ci£, 1 + C2£ 



!+/£,! + Ci£, 1 + C2£ 

The two bases defined by Eqs. (122) and (123) are related as 



(123) 



4-^3 I 1 



ai£, a2£, a3£, a4£ 

2 + /£, 1 + Ci£, 1 + C2£ 



1 + 



2aia2a3a4£ 



(l + 2/£)(l + Ci£)(l + C2£) 



5i^4 



^4i^3 



= 2z 



ai£, a2£, a3£, a4£ 

i + /£, 1 + Ci£, 1 + C2£ 
4 

a 1 020.3 0.4 £ 

(l + 2/£)(l + Ci£)(l + C2£) 



4-^3 



1 + ai£, 1 + a2£, 1 + a3£, 1 + a4£, 1 

I + /£, 2 + Ci£, 2 + C2£, 2 



1 + ai£, 1 + a2£, 1 + a3£, 1 + a4£ 

I +/£,2 + Ci£,2 + C2£ 



.(124) 



The second derivative in Eq. (122) may be obtained from any of the following 
relations: 



4i^3 

4-P3 



1 + ai, 1 + 02, 1 + 03, 1 + 04 

1 + /, 1 + Ci,C2 

1 + Oi, 1 + 02, 1 + 03, 2 + 04 
1 + /,1 + Ci,l + C2 



fcie{e + c2-i) /ai, 02, 03, 04 
4^3 f 

Z 01020304 \ J,Cl,C2 

/ ciC2e{e + 04) ^ /ai, 02, 03, 04 

201020304(1 + 04) ^\ f,Cl,C2 



z 

(125) 
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with the substitutions of Eq. (116). In particular, it can be expressed it as a hnear 
combination of the second, third, and fourth hypergeometric functions in Eq. (123) 
as 



(71+72) /C1C2 

z 01020304 

- [0472 + (C2 - 1)71] 4-P3 ^ 
+72(1 + 04) 4F3 



71 C2 4i^3 



1 + Oi, 1 + 02, 1 + 03, 1 + 04 
1 + /, l + Ci,C2 



Ol, O2, 03, 04 
/, Cl,C2 

1 + Oi, 1 + 02, 1 + 03, 1 + 04 
1 + /,1 + Ci,l+C2 

1 + oi, 1 + 02, 1 + 03, 2 + 04 

1 + /,1 + Cl,l + C2 

where 71 and 72 are arbitrary numbers. Putting 72 = —71, we obtain a hnear 
relation between the hypergeometric functions on the r.h.s. of Eq. (126). We checked 
that the e expansions of the hypergeometric functions constructed in Ref. [20] satisfy 
this relation. The third derivative in Eq. (122) may be obtained from any of the 
following relations: 



(126) 



4-^3 

4i^3 

4-^3 



1 + Oi, 1 + 02, 1 + 03, 1 + 04 
1 + /,Cl,C2 

1 + oi, 1 + 02, 2 + 03, 2 + 04 

1 + /,1 + Cl,l + C2 

1 + Oi, 1 + 02, 1 + 03, 2 + 04 

1 + /, 1 + Ci,C2 



/^(^ + Ci-1)(0 + C2-1 



/ CiC2g(g + 03)(g + 04) 
Z 01020304(1 + 03)(1 + 04) 
/Cig(g + 04)(g + C2-l 

z 01020304(1 + 04) 



4-^3 
4i^3 
4i^3 



«!, «2, 03) «4 
/, Cl,C2 

01)02)03) O4 
/) Cl,C2 

Ol) 02)03)04 
/, Cl)C2 



(127) 
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with the substitutions of Eq. (116). As a consequence of Eq. (127), the fohowing 
Unear relations between hypergeometric functions are vaUd: 



Cl 4-^3 

= (Ci - 04 - 1) 4F3 
C2 4-P3 

= (C2 - 03 - 1) 4F3 
(/3l+/32)<CiC2 4i^3 



(1 + 04) 4i^3 



1 + Oi, 1 + 02, 1 + 03, 1 + 04 

1 + /, Cl,C2 

1 + Oi, 1 + 02, 1 + 03, 1 + 04 
1 + /, 1 + Ci,C2 

1 + oi, 1 + 02, 1 + 03, 2 + 04 

1 + /, 1 + Cl,C2 

1 + oi, 1 + 02, 1 + 03, 2 + 04 

1 + /,1 + C1,1 + C2 



(1 + 03) 4i^3 



1 + Oi, 1 + 02, 1 + 03,2 + 0,4 
1 + /, 1 + Ci,C2 



1 + oi, 1 + 02, 2 + 03,2 + 04 

1 + /,1 + Cl,l+C2 



1 + oi, 1 + 02, 1 + 03, 1 + 04 

1 + /, Cl,C2 



1 + oi, 1 + 02, 2 + 03, 2 + 04 

1 + /, 1 + Cl, 1 + C2 



-(l + 03)(l + 04) 4i^3 

+ (C2 - 04 - 1) [Pl{c2 - 03 - 1) - P2{C1 - 04 - 1)] 4F3 
-,01 (Cl + C2 - 2 - 03 - 04)02 4-^3 
-/02(C1 + C2 - 2 - 03 - 04)(1 + 04) 4F3 



1 + Oi, 1 + O2, 1 + O3, 1 + 04 
1 + /,1 + Ci,l+C2 

1 + Oi, 1 + 02, 1 + 03, 1 + 04 

1 + /, 1 + Cl,C2 

1 + oi, 1 + 02, 1 + 03,2 + 04 

1 + /,1 + Cl,l + C2 



0, 



(128) 



where Pi and P2 are arbitrary numbers. We checked that the e expansions of the 
hypergeometric functions 4F3 constructed in Ref. [20] and collected in Ref. [48] 
satisfy all these relations. 

In the case oi = 1, the basis of the differential reduction is 



1,1 + 0,26, 1 + 036, 1 + 046 



f + /e,2 + cie,2 + C2e 
The first derivative in Eq. (129) can be related to any of the two functions. 



(129) 



4i^3 



1,1 + 026, 1 + 036, 2 + 046 

I + /£, 2 + Ci£, 2 + C2£ 



4-^3 



1,1 + 026, 1 + 03^, 1 + 046 
I + /£, 2 + Ci£, 1 + C2£ 



(130) 

by single application of the operators B^^^ or Hc^ ^ given in Eq. (7). The second 
derivative in Eq. (129) is related to any of the following hypergeometric functions: 



4-P3 

4i^3 



1 + ai£, 1 + 02£, 1 + a3£, 1 + 04£ 
I + /£, 1 + Ci£, 1 + C2£ 



4i^3 



1 + ai£, 1 + 02£, 2 + 03£, 2 + 046 
I + /£, 2 + Ci£, 2 + C2£ 



1 + Oi£, 1 + 02£, 1 + 03£, 2 + 04£ 
3 
2 



I + /£, 2 + Ci£, 1 + C2£ 



(131) 



which are generated from the basis function of Eq. (129) via application of H^^ ^i^^^ \ 



-Das ^04 ' -"63 • 
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A.2.4 5F4 

Similarly, for the hypergeometric function 



h + aie, I2 + a2e, h + ase, I/^ + a^e, Is + 056 
-^6 + ^ + /e, h + cie, /g + C2e, /g + £36 



(132) 



where {1^} are arbitrary integers and all 7^ 0, the basis of the differential reduction 



IS 



aiE, 026, 036, a4e, ose 



5 + /e, 1 + cie, 1 + C2e, 1 + 036 



(133) 



In this case, there are eleven functions of the type of Eq. (109). The two bases 
defined by Eqs. (133) and (109) are related as 



04s, a.2£, Cl3£, (346, CL^£ 

i + /e, 1 + cie, 1 + C2e, 1 + c^e 
2e^ 



z]=l + 



(n ■=! a) 



(1 + 2/e) nLi(l + Cik£) 

^ 5-^4 



1 + aie, 1 + 026, 1 + 035, 1 + 046, 1 + 056, 1 

I + /£, 2 + Ci£, 2 + C2£, 2 + C3£, 2 



ai£, a2£, 036, 046, a^e 

i + /£, 1 + Cl£, 1 + C2£, 1 + C3£ 



2z 



(nLi«i) 



5i^4 



(i+2/£) [nLi(i+cfc£ 

The second derivative in Eq. (133) may be obtained as 

(71 + 72) ^ (nLi 



1 + ai£, 1 + a2£, 1 + a3£, 1 + a4£, 1 + a^e 

I + /£, 2 + Ci£, 2 + C2, 2 + C3££ 



(134) 



0,1, 02, 03, 04, as 

/, Cl,C2,C3 



= 71C3 5-P4 

+ 72(1 + 05) 5i^4 
- [71 (C3 - 1) +72^5] 5-F4 



1 + ai, 1 + 02, 1 + 03, 1 + 04, 1 + as 

1 + /, 1 + Ci, 1 + C2,C3 

1 + ai, 1 + 02, 1 + 03, 1 + 04, 2 + 05 

1 + /£, 1 + Ci, 1 + C2, 1 + C3 

1 + Oi, 1 + 02, 1 + 03, 1 + 04, 1 + 05 
1 + /£, 1 + Cl, 1 + C2, 1 + C3 



, (135) 



with the substitutions of Eq. (116), where 71 and 72 are arbitrary numbers. Putting 

72 = —71, we obtain a linear relation between the hypergeometric functions on the 
r.h.s. of Eq. (135). We checked that the e expansions of the hypergeometric functions 
constructed in Ref. [20] satisfy this relation. The third derivative in Eq. (133) is 
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related to one of the following functions or their linear combination: 



5-^4 

5-^4 
5-^4 

where 



1 + ai, 1 + a2, 1 + 03, 2 + 04, 2 + 05 

1 + /,1 + Ci,l + C2,l + C3 

1 + ai, 1 + 02, 1 + as, 1 + 04, 2 + 05 

1 + /, 1 + Ci, 1 + C2,C3 

1 + ai, 1 + 02, 1 + 0,3, 1 + a4, 1 + 05 

1 + /, l + Ci,C2,C3 



f[Ut=iCk) eie + a,)ie + a5) 

' (n|=i«i)(l + «4)(l + a5) 
/(nLiCik)^(^ + C3-l)(^ + a5) 



n-=i% (1+05) 



f Ci9ie + C2 - l)i9 + C3 - 1) 



(n5=ia.) 



F{z)= 5F4 



Ol, 02, 03, 04, 05 
/, Ci,C2,C3 



, (136) 



(137) 



As a consequence of Eq. (136), the following linear relations between hyper geometric 
functions hold: 



(1 + 04 - C3) 5F4 
(1 + 04) 5F4 

- C3 5^4 

(1 + as - C2) 5-P4 
(1 + as) 5^4 

- C2 S-^4 



1 + ai, 1 + a2, 1 + as, 1 + 04, 2 + as 

1 + /,1 + Ci,l + C2,l + C3 

1 + ai, 1 + a2, 1 + as, 2 + a4, 2 + as 

1 + /,1 + Ci,l + C2,l + C3 

1 + ai, 1 + a2, 1 + as, 1 + 04, 2 + as 

1 + /, 1 + Ci, 1 + C2,C3 

1 + ai, 1 + a2, 1 + as, 1 + 04, 1 + 05 

1 + /, 1 + Ci, 1 + C2,C3 

1 + ai, 1 + a2, 1 + as, 1 + a4, 2 + as 

1 + /, 1 + Ci, 1 + C2,C3 

1 + ai, 1 + a2, 1 + as, 1 + a4, 1 + as 

1 + /, l + Ci,C2,C3 

(as - cs + 1) [Pi (cs - 04 - 1) - /32 (c2 - as - 1)] 5F4 

= {Pi + /32)C2C3 si^4 

-(^i+^2)(l + a4)(l + as) si^4 
+ P1C3 (a4 + as - C2 - Cs + 2) s^4 
+ /32(1 + as) (a4 + as - C2 - C3 + 2) 5F4 



1 + ai, 1 + a2, 1 + as, 1 + a4, 1 + as 

1 + /,1+Ci,l + C2,l + Cs 



1 + ai, 1 + a2, 1 + as, 1 + 04, 1 + as 

1 + /, l + Ci,C2,C3 

1 + ai, 1 + a2, 1 + as, 2 + a4, 2 + as 

1 + /,1 + Ci,l+C2,l+Cs 

1 + ai, 1 + 0,2, 1 + as, 1 + a4, 1 + as 

1 + /, l + Ci,l + C2,Cs 

1 + ai, 1 + a2, 1 + as, 1 + 0,4, 2 + 05 

1 + /,1 + Ci,l+C2,l + Cs 
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The forth derivative in Eq. (133) is related to any of the following hypergeometric 
functions: 



5-^4 
5F4 



1 + ai,l + 02,2 + 03,2 + 04,2 + 05 

1 + /,1 + Ci,l + C2,l + C3 

1 + Oi, 1 + 02, 1 + 03, 2 + 04, 2 + 05 

1 + /, 1 + Cl, 1 + C2,C3 

1 + Oi, 1 + 02, 1 + 03, 1 + 04, 2 + 05 

1 + /, l + Ci,C2,C3 

1 + Oi, 1 + 02, 1 + 03, 1 + 04, 1 + 05 
1 + /, Ci,C2,C3 



/ 0102036 {6 + as) je + a4){d + ab) 
'(UU^'j) (l + a3)(l + a4)(l + a5) 

/ciC2g(g + C3-l)(g + 04)(g + Q5) 

' (n5=ia,)(l + a4)(l + 05) 

f ci9{e + 02 - i){e + 03 - i){e + a5 



F{z) 



n •=! «i j (1 + 05) 

f 0(6 + 01 - i)i0 + 02- i){e + 03 -I) 

z 



(n5=i«.) 



Fiz), 
(139) 



with the substitutions of Eq. (116). The set of algebraic relations between the 
hypergeometric functions of Eq. (139) is to long to be presented here. We checked 
that the e expansions of the hypergeometric functions 5F4 constructed in Ref. [20] 
satisfy all these linear relations. 



A.2.5 p+iFp, p>6 

For the hypergeometric function 



(140) 



the following representation was derived in Refs. [67]: 



- P 



2'{2}p-l 



^ {2n 

k 



1 



1 in) ^ 

p-2 



i=o 



{-2y 



[a 



{Inzf-^-^Lsfl^ { 2arcsin^ ) (141) 



where < ^ < 4, p > 2, and Ls^ (z) denotes the generalized log-sine function^^ [69]. 
Later, it was shown in Ref. [29] that any generalized log-sine function of argument 
9 can be written in terms of Nielsen polylogarithms with argument^^ y = exp(±i^) 
(see Eqs. (2.18) and (A.21) in Ref. [29]). 



^^A computer program for the numerical evaluation of the generalized log-sine functions with high 
precision was developed in Ref. [68]. 

^^In terms of the variable z in Eq. (141), we have 



z /— z 
V = COS0 ± isin^L „ . ^=l — -±iy/zxl — 7 = 

» le=2arcsin^ 2 V 4 ]^ ^ 
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